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SUMMARY 


The application of a combined linear theory/impact theory 
method to calculate pressures and loadings on a wing^body 
configuration at Mach 4.63 was assessed. The results, compared 
with experimental pressure data, show that the combined method 
gives significantly improved predictions over either the linear 
theory or the impact theory method alone. The combined method was 
also applied in the inverse design mode at Mach 4.63 to calculate 
optimum camber distributions on a wing-alone, and on a wing-body 
configuration. The results of these optimization calculations are 
compared to results obtained using unmodified linear theory, and 
show that there is a large difference in the predicted camber 
distributions. Assessment of the analysis results indicates that, 
for the high Mach number, the optimum camber slopes obtained using 
the combined method are more correct than the linear theory 
results, and that finite thickness wings "optimized" at high Mach 
numbers using unmodified linear theory will not achieve the 
minimum drag characteristics for which they are designed. 


INTRODUCTION 


The feasibility of combining elements of linear theory and 
impact theory for improved aerodynamic predictions in the Mach 4 
to 8 range was verified in Reference 1 for several wing-alone 
configurations. The approach is to use the aero influence coeffi- 
cients (AIC) from the linear theory computer program of Reference 
2 to modify the impact pressures computed in the impact theory 
analysis program of Reference 3. The approach provides a method 
of accounting for the interference effects missing from the impact 
theory analysis, while retaining the thickness and non-linear 
lifting effects characteristic of the high Mach number range. 



The basic equation for the combined theory, given in matrix 
form, is 


-1 g * 

{Cp} = [a 1 { 4 Cp} 


( 1 ) 


where i Cp } is the column matrix of upper or lower surface pres- 
sure coefficients, [a“l] is the square matrix of aero influence 
coefficients, g = /m^- 1 where M is the freestream Mach number, and 
Cp* represents the upper or lower surface pressure coefficients 
calculated using impact theory. The approach requires a one-to- 
one correspondence between the linear theory chord plane panels 
used to generate the AIC matrix and the surface elements used to 
represent the configuration for the impact theory calculation, 
i.e., the projection of the surface elements onto the chord plane 
must match the linear theory panels. As noted in Reference 1, the 
impact pressure options consistent with this analysis are tangent 
wedge, and Prandtl-Meyer . 

The combined theory was applied to several wing-alone config- 
urations in Reference 1, with excellent results for pressures, 
loadings, and forces and moments. Only in the case of comparisons 
near the outboard leading edge of a 76 degree swept wing were 
there substantial differences between predicted and measured 
pressures. Further analysis of that problem indicates that the 
differences are related to the high sweep angle, and that the 
over-predictions encountered will not occur unless the wing sv/eep 
exceeds about 70 degrees. 

In Reference 1, the basic equation (1) was used as the start- 
ing point for a derivation of closed form approximations for the 
lift, drag, and moment on uncambered, wing-like bodies at high 
Mach numbers. Comparisons made with force and moment data on 
several NASA all-body models, at Mach 5.37 and 7.38, show that 
these closed form equations correctly predict the manner in which 
thickness, volume, and non-linear angle of attack effects modify 
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the linear theory predictions. Because of these effects, 
unmodified linear theory applied in the inverse optimization mode 
will give camber slopes which correspond to incorrect lift and 
moment constraints. In Reference 1, the combined theory was used 
to rederive the optimization equations to properly account for the 
thickness and higher order effects on the camber distributions. 

The application of the combined theory to a wing-body config- 
uration to calculate pressures and loadings, and to calculate 
optimum camber distributions, is discussed in this report. The 
necessary modifications to allow the computer programs given in 
References 2 and 3 to be used for the combined theory are also 
discussed . 


WING-BODY ANALYSIS 


The wing-body configuration discussed in Reference 4, and 
shown in Figure 1, was analyzed at Mach 4.63 using the combined 
analysis. The 65 degree swept delta wing has a symmetrical double- 
wedge airfoil section of 6 percent thickness ratio and joins the 
slender axisymmetric body at 47 percent of the body length. 

Pressure data is given on the body and at four spanwise locations 
on the wing. 


For a wing-body, the basic equation (1) for the combined 
analysis can be separated into two equations, one for the body and 
one for the wing. To do this, the AIC matrix (a“^) is subdivided 
into the submatrices which represent the body-on-body panel 

influence (agl), the wing-on-body influence (a^l), the body-on- 

wing effects and the wing-on-wing panel effects (a^). The 

basic equation then becomes 



d') 
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In terms of lifting pressures, the equations for the body and wing 
are, respectively, 

( 2 ) 

= '^BwJ 4 I 

In the equation for the body pressures (^Cp^) , the first term is 

simply the isolated, or body alone, pressure, and the second gives 
the interference contribution from the wing. In the equation for 
wing pressures (ACp^), the first term is the interference contribu- 
tion from the body, and the second is the wing alone solution. 



Figure 1. Sketch of Wing-Body Modei from NASA TN D-6480. 
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Two linear theory methods were considered for the calculation 
of the influence coefficients. One was the program described in 
Reference 5, which uses surface panels on the body and chord plane 
panels for the wing. Although this program gave excellent predic- 
tions for the body alone pressures, a stable solution in the 
presence of the wing could not be achieved. 

The second approach was to use an all chord plane represen- 
tation (i.e., flat panels of zero thickness), as illustrated in 
Figure 2, and to use the program of Reference 2 to compute the 
influence coefficients. The pressures and loadings on the body 
were calculated for a variety of impact pressure options using 
these chord plane influence coefficients, and the results, partic- 
ularly on the forebody, v;ere poor. Figure 3 presents typical 
loadings on the body calculated using the chord plane influence 
coefficients. For comparison, the surface panel solution and the 
experimental pressure differences along the centerline are also 
shown. The large pressure peaks shown in Figure 3 are attributed 
to the high sweep (81 to 86 degrees) of the leading edge panels. 

On the afterbody, the loadings become more reasonable. 
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Figure 3. Predicted Loadings on Body Using Chord Plane Panels 

Mach 4.63 Angle-of-Attack = 6.2 Deg 






For the body pressure calculation, it is not really necessary 
to use the ag^ influence coefficients as computed. If the body- 

on-body influence coefficient matrix is replaced by a diagonal 
matrix, where the diagonal terms are 4/g, the combined theory will 
give the impact theory solution on the body, plus the wing carry- 
over effects. If the [^ggl matrix is zeroed out entirely, then 

only the wing interference effects are retained. These can be 
calculated and added to a separate isolated body solution, e.g., 
the surface panel method of Reference 5. 

Figures 4 and 5 show the predicted pressures and loadings on 
the body using the latter approach. The isolated body pressures 
were computed using the surface panel method of Reference 5, and 
the wing interference pressures, calculated using tangent wedge 
and Prandtl-Meyer pressure options, were added directly. These 
results, compared with the experimental data, show that both the 
pressure levels and the loadings are very well predicted, both on 
the forebody, and in the interference region on the afterbody. 

The results of the chord plane body representation shown in 
Figure 3 indicate that the forebody effects on the wing will be 
overpredicted. On the other hand, since the results on the after- 
body were more reasonable, the carry-over effects from the panels 
aft of the wing-body junction should also be reasonable. 

The pressure distributions on the wing presented in Figures 6 
through 9 were computed using tangent wedge/Prandtl-Meyer (TW/PM) 
impact pressures on the wing and aft portion of the body. Three 
different pressure options were considered on the forebody. These 
were: (1) tangent wedge/Prandtl-Meyer, where the forebody was 

treated as a flat plate. This calculation gives essentially the 
same carry-over loadings on the wing that unmodified linear theory 
gives. (2) A tangent cone (TC) approximation was used, where 

C* = + 2 sin^6 
^TC 
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x/e 


Figure 4. Pressure Distribution and Loading on Body Centerline, 
Mach 4.63 Angle-of-Attack = 6.2 Deg 



— Isolated body solution 

Isolated body + wing 
carry-over contribution 

O Data (NASA TN D-6480) 



Figure 5. Pressure Distribution and Loading on Body (60 and 120 Deg Meridians), 

Mach 4.63 Angle-of-Attack = 6.2 Deg 






Figure 6. Wing Pressures and Loadings (y/b/2 = 0.258), 
Mach 4.63 Angle-of-Attack = 6.2 Deg 
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Figure 7. Wing Pressures and Loadings (y/b/2 = 0.444), 
Mach 4.63 Angle-of-Attack = 6.2 deg 
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and 6 is the local surface slope measured from the freestream. 

For negative values of s, the minus sign in the above equation is 
used. And (3), the wing pressures were computed with the forebody 
pressures set to zero. 

The results of the calculations for the wing pressures are 
compared with the measured data in Figures 6 through 9. Results 
from unmodified linear theory on a chord plane representation of 
the wing are also shown, along with a standard impact theory 
(TW/PM) solution for the wing surface panels. The comparisons at 
all four wing stations show that the combined theory gives better 
predictions for both pressures and loadings than either linear 
theory, or impact theory. For the combined theory results, the 
best predictions were obtained for the forebody pressures set to 
zero, suggesting that the slender forebody actually contributes 
very little influence on the wing. The use of both the tangent 
cone approximation, and the tangent wedge/Prandtl-Meyer 
calculation, led to overpredictions at the inboard wing stations. 
These overpredictions occur along the Mach line emanating from the 
highest pressure peaks on the forebody (Figure 3), and reflect the 
difficulty in simulating the flow on a slender axisymmetric body 
using a chord plane representation. 


OPTIMIZATION 

Figures 6 through 9 show that at each spanwise location on 
the wing, the center of lift is much further forward than 
predicted by linear theory. This is because the thickness contri- 
butions to the lifting pressures, which become important at high 
supersonic Mach numbers, are not included in the linear theory. 
Depending on the thickness distributions, the linear theory will 
not correctly predict the lift and/or moment characteristics of 
finite thickness wings at high Mach numbers. Conversely, the 
"optimum" camber distributions calculated using unmodified linear 
theory will correspond to incorrect lift and moment constraints. 
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As discussed in Reference 1, the combined theory permits the 
inclusion of the thickness effects in the inverse design proce- 
dure. The optimization equations from Reference 2 were modified 
to include the thickness contributions to the lifting pressures, 
and the equations are presented in Figure 10. These equations 
were used to calculate the optimum camber slopes on the theoreti- 
cal wing from the wing-body configuration in Figure 1. The camber 
slopes were calculated at Mach 4.63 for a design lift coefficient 
(Cl) of 0.1 and a design moment coefficient (0^) of 0.0. The 
moment reference center was at 50 percent of the mean aerodynamic 
chord (.5c). The calculated camber slopes using the combined 
theory are compared with the unmodified linear theory results in 
Figure 11. 

The comparisons in Figure 11 show that there is a significant 
difference in the predicted camber distributions. The reason for 
the large difference is illustrated in Figure 12, which shows the 
lift and moment characteristics for the uncambered wing as 
predicted by the combined theory and linear theory. Although the 
lift curve slope for this wing is the same for both theories, the 
combined theory shows that the wing is 9.5 percent unstable, while 
the linear theory predicts the wing to be neutrally stable. Thus, 
the moment constraint (Cji^ - = 0.0) is not a driver in the linear 
theory analysis, and the camber slopes do not deviate much from a 
6.63 degree flat plate angle of attack needed to obtain a lift 
coefficient of 0.1. 

On the other hand, the combined theory predicts that, at a 
lift coefficient of 0.1, the camber must effectively trim out a 
moment coefficient of 0.0095. Since the thickness effects shift 
the loading forward, relative to the linear theory, the forward 
portion of the wing must be drooped more than predicted by the 
linear theory to decrease the moment, and to satisfy the moment 
constraint. As a consequence of the greater variation in camber 
slopes, the minimum drag due to lift is about 5 percent higher 
than predicted by the linear theory alone. 
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OPTIMIZATION MATRIX - Solution gives loadings for minimum drag 
corresponding to lift and moment constraints 
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OPTIMUM CAMBER - Calculated from optimum loadings 



1 s 

“i = ~ j ^ij <^Cp.)opT 

where Sj 

= area of element i 

"ij 

= influence coefficient (element i on element j) 

bi 

= 1+1 .2/3 Cj + 0.6)3^^^ €j^ (thickness factor from 
Ref. 3 Program) and 6j is thickness slope 


= pressure difference on element i 

X — Xj 

= distance from moment ref to centroid of element 

Ai, A2 

= Langrangian multipliers 

r, M 

= design lift and moment constraints 
L = C[_ qS M = C|^ qSc 


Figure 10. Optimization Matrix for Design Procedure. 



“camber deg “camber deg 


V/b/2 = 0.05 


V/b/2 = 0.25 



x/c 



0 0.2 0.4 0.6 0.8 1.0 

x/c 


— Combined theory 
(Cq = 0.0115) 

— — — Unmodified linear 

theory (Cq = 0.01 1 0) 


6% Double-Wedge Airfoil 




x/c 




0 0.2 0.4 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0 

x/c x/c 


Figure 11. Predicted Camber Distributions for Minimum Drag 

65“ Swept Delta Wing @ Mach 4.63 
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Figure 12. Lift and Moment Characteristics of Uncambered Wing 
65“ Swept Delta Wing at Mach 4.63 
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The minimum drag characteristics for the wing-body (Figure 
1) , for the same lift and moment constraints used above, are illus- 
trated in Figure 13 as a function of the moment reference location 
(Xcg)* these calculations, the panels representing the body 

were constrained at 6.2 degrees angle of attack. Without the 
camber constraint on the body, the solution tends to give unreal- 
istically large camber slopes on the forebody. 
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Figure 13. Minimum Drag Due to Lift as Function of Moment Reference Location. 

Wing-Body Configuration @ Mach 4.63 
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Assuming the combined theory results reflect the true situa- 
tion (a reasonable assumption in view of the results shown in 
Figures 6 through 9), the minimum drag is achieved for the center 
of gravity at 34 percent of the mean aerodynamic chord. For this 
moment reference location, the linear theory slopes give higher 
drag. Since the neutral point is incorrectly predicted, the 
linear theory camber slopes will also fail to meet the moment 
constraint. Because of this, trimming devices would have to be 
used, resulting in even higher drag. 


COMPUTER PROGRAMS 


The computer programs described in References 2 and 3 were 
used to perform the combined linear theory/impact theory calcula- 
tions described in this report and in Reference 1. Although 
various versions of the Reference 2 program may be in use, so long 
as the essential subroutines (EVAL, INVW, MDMATE, FCALC, CAMBW and 
DCPD) are contained in those programs, they may be used for the 
linear theory portions of the combined calculation. For these cal 
culations, only minor modifications to subroutines INVW, MDMATE, 
CAMBW, and DCPD are required. It should be noted that only INVW 
need be modified to obtain the inverted AIC matrix for use in the 
impact theory analysis. For optimum camber calculations, MDMATE, 
CAMBW, and DCPD need to be modified. The actual modifications to 
these routines are given in Appendix A. 

The impact theory program (Reference 3) is used to calculate 
the impact pressures, and, using the AIC matrix, to calculate the 
combined theory pressures. To perform the pressure calculation, a 
new subroutine (CPCALC) is added to the program, along with minor 
changes to subroutines PRES and FORCE which provide required 
control information. The thickness factors needed in the optimize 
tion matrix are calculated in CPCALC. The modifications to the 
program are given in Appendix B. 
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RESULTS AND CONCLUSIONS 


The results of this study demonstrate that the combined 
linear theory/ impact theory provides an accurate, versatile, and 
easy to use method for the analysis and design of configurations 
at high Mach numbers. The combination of the supersonic linear 
theory and the hypersonic impact theory overcomes many of the 
problems encountered by the individual theories when used in the 
Mach 4 to 8 range. In addition, only the simplest combination of 
characteristics from the linear and impact approaches are required 
to produce consistently good results. Oblique shock theory (or 
Newtonian) and Prandtl-Meyer expansion coupled with a chord plane 
representation of the geometry to generate the aerodynamic 
influence coefficient matrix is considered the basic calculation 
method. However, by partitioning the aerodynamic influence coef- 
ficient matrix, the basic approach becomes flexible enough to 
permit a multitude of combinations which may become more evident 
through continued applications and comparisons with experiment. 

The results of the wing-body comparisons show that the 
combined theory gives improved analysis results over either the 
linear theory or the impact theory methods alone. The application 
of the combined theory in the inverse design mode shows that the 
use of unmodified linear theory to calculate optimum camber slopes 
will lead to incorrect results at high supersonic Mach numbers. 

The combined theory offers the prospect of improved high speed 
designs by providing a more appropriate optimization procedure. 
Only minor modifications to existing state-of-the-art computer 
programs are needed to implement the combined theory, and these 
modifications will not affect the normal, stand-alone operation of 
the programs, nor require changes in existing input data decks. 


20 



APPENDIX A 

LINEAR THEORY PROGRAM MODIFICATIONS 

The modifications to the linear theory program documented in 
NASA CR-73107 (Part II) are given below. The listings of the sub- 
routines affected have been extracted from the NASA documentation 
and the appropriate changes are noted. 

In the analysis mode, the program is used to generate the AIC 
matrix and to write the matrix onto a tape, which can be accessed 
by the impact theory program. It is assumed that the user will 
specify an appropriate output device and the necessary job control 
cards to save the output tape. Further, it is assumed that the 
AIC matrix will be calculated using the "wing alone" or all chord 
plane paneling option. The AIC matrix is obtained from subroutine 
INVW, as shown in Figure A-1. Note that this is the inverted 
matrix, [a“^] . 

For the inverse design or optimum camber calculations, the 
MDMATE and CAMBW routines are modified to read in and to use the 
thickness factor terms from the impact theory program. Again, it 
is assumed that the user will supply the appropriate tape number 
and job control cards to allow the program to access the output 
tape from the impact theory program. Subroutine DCPD should also 
be modified in case the configuration is run through a series of 
angle of attacks using the calculated camber slopes. The 
modification in DCPD corrects the loading distributions (ACp's) 
consistent with the combined theory and gives the corrected lift 
and moment characteristics. The calculated camber distributions, 
and the corresponding loadings are printed in the normal output 
(pages 102-104 in Part II) under the headings "WING CAMBER SLOPES 
(DZ/DX)" and "WING PANEL PRESSURE DIFFERENCE (CL)," respectively. 
Note that with the exception of the input of an indicator to tell 
the program to read a tape, all other inputs are unchanged. If 
the indicator is placed in an unused field on one of the standard 
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SUBROUTINE INVH 


C FOR WING ONLY CASE. STORES MATRIX (A) AND tAt INVERSE ON TARE 

C MAXIMUM SIZE MATRIX INVERSION >110 


COMMON DATE! Z) .NTAPEA.NTAREB.NTAREC.NTARED.NTAREE.NTAREF.NTAREI 
l.NTAPEO.NBODY.NWING.XMACH.SVM.KACE 

DIMENSION AUHdlS.ltSI 

MCEMEN-115 

C READ <A) MATRIX INTO CORE. WRITE ON TARE 

CALL FSF( 1.NTAPEA.IRR1 
00 100 J'l.NUING 

READ (NTAPEA) { AWWI I.J).I«1,NWINGI 
WRITE INTAPEEi I AWW( I . J 1. 1 >1.NWING» 

100 CONTINUE 

ENO file NTAPEE 
REWIND NTAPEA 


C INVERT MATRIX (Al 

CALL SINVRTI AWW.HOEHEN.NWINC. IRR1.IRR2, SCALE. OET.NOETXRI 
IF (IRRII ISO.ZOO.IBO 
150 WRITE INTAPE0.6000) IRRl. IRRZ. SC ALE 
REWIND NTAPEE 
STOP 


FROM PAGE 671 OF 
REF. 2 - PART II 



200 CONTINUE 

DO 750 J«I.NW1NG 

^ WRITE (NTAPEE) ( AWMI I . J I. I >1.NWINCI 
250 continue 

ENO file NTAPEE 
REWIND NTAPEE 


*000 FORMAT! 1H1,38HERB0R IN INVERSION OF MINC ONLY MATRIX 
1.5X.0HIRRI >. 13. 5X.6HIRR2 >.I3.5X,7HSCALE >.E12.6) 

RETURN 

ENO 



MTRPE£= ■ » 

WRITE <NThPES) I> > 1 = 1 » Hl.lIWG) 


* REQUIRES DEFINITION OF NTAPES 


Figure A-1. Modifications to Subroutine INVW. 
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input cards r such that it defaults to zero if not specified, then 
the modifications, given in Figures A-2 to A-4, will not affect 
the normal operation of the program. 


SUBROUTINE NDB4TE 



C FORM»TION OF DRAG NININIZATION NATRIX 

C XAXIMUM SIZE MATRIX INVERSION -IIZ 

C MUST BE 2 PLUS MAXIMUM NUMBER OF WING PANELSIIIOI 

COMMON OAT El ZI.NTAPEA.NTAPEB.NTAOEC.NTAPEO.NTAPEE.NTAPEF.NTAREI 
l.NTAPEO.NBOOY.NWlNG.XMACH.SYn.KACE 


DIMENSION WWI IIS, IISI.XBARI ZIOi.AREAl 2101 
M0FMFN»115 

READ (NTAPECI NBOOY, NWING. XMACH, SYM, RACE 
NPANEL-NBOBY*NWING 

READ (NTAPECI I I.XBAR(II.0UHMY1,DUMMY2,DUMHY3,DUNHY«,DUHMYS 
I < AREA! I I .DUMMY 6f OUMMYT, OUMMYA, I » l.NP ANEL I 
REWIND NT APEC 



DO too J*1,NWING 

READ (NTAPEEI I WWI I , J I. I< I .nWINGI 

100 CONTINUE 

REWIND NTAPEE 


DO 200 I>I,NWING 
I-I 

ii»i»nbooy 

DO 200 J=l, I 
JJ- J«NBODY 

WWI I ,J1=-WW( I, JIPAREAf II)-WW(J,II*AREA( JJI 
WWIJ.I I.WWI I.JI 
200 CONTINUE 


FROM PAGE 686 
REF. 2 - PART 


OF 

II 



C: □ M f'1 □ N E; F H C T ••• B M D ■ 1 1 0 :• 



DD jij J=1 n NMIHb 
EMD '■ j;:' = 1 . U 


PERU ':;5' 99 ' BIND 
99 FDPMRT ':'F1 0. 1> 

IF '::BINri. EG!. 0. O;:' GO TO l?n 
HTRPET= 

PERU 'rMTHPET ' •TBMDGlJ i. I = 1 1. H!.l I HG :> 
DO 15 0 .1=1 !. r-(l..lIHG 
DO 150 I = 1 Hl.l I NG 
1 5 0 i.illii 1 1 _l :> =liH,il C I . J ;■ ■•• Er'ID I :> 

170 CONTINUE 


* REQUIRES DEFINITION OF NTAPET 


Figure A-2. Modifications to Subroutine MDMATE. 
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0 


SUBROUTINE CAMBW ( NW .NT APEX , A.CLW.ALPHAW ) 
DIMENSION A(l'),ALPHAwm,CLh(ll 


FOR MING ONLY C ASE COMPOTE S NORMAL VELOCITY COMPONENTS ON WING 


100 



CO ICO J=1.NM 
ALPHAWI Jl=0.0 
CONTINUE 
DO 200 J-l.NW 

READ (NTAPEX) (A( I I ,1 = 1.NU) 

00 200 1=1. NW 

ALPHAWI I )=ALPHAW( I I+AI I )*CLU1J) 
CONTINUE 

RETURN 

END 


FROM PAGE 506 OF 
REF, 2 - PART II 



COMMON .•' EFRCT -'EMD O 1 0) 



DO 3i;ii:i 1 = 1 jNw 
RLPHRI.I C I ) =RLPHRI.) < I > /BMD -rn- 
Snn CONTINUE 


Figure A-3. Modifications to Subroutine CAMBW. 



SUBROUTINE DCPDI NM, NTAPEX , A, ALPHAH, CIH) 
COMPUTES wing PANEL PRESSURE DIFFERENCE 
DIMENSION Al 1 ) , ALPHAMtl ) tCLMI 1) 



DO ICO J=1,NM 
CLMI J)=0.0 
100 CONTINUE 

DO 200 J=1,NM 

READ (NTAPEX; (Ail), 1=1, NM) 
DO 200 1=1, NH 

200 continue 

return 

END 


FROM PAGE 535 OF 
REF. 2 - PART II 



COMMON /E;FRr:T.''EMD Q 1 0> 



CLMO'; 


=CLM a") +R a:> ♦rlphrm cj;.. ♦bmd 




Figure A-4. Modifications to Subroutine DCPD 



APPENDIX B 

IMPACT THEORY PROGRAM MODIFICATION 


The modifications to the impact theory program documented in 
AFFDL-TR-73-159 (VOl. Ill) are discussed below. Two subroutines 
(PRES and FORCE) are modified, and a new subroutine (CPCALC) has 
been added. 

In the PRES routine, two control parameters (IWOOD and IROW) 
are read in and placed in common. The required changes are given 
in Figure B-1. The card identification numbers correspond to the 
ID numbers in the reference listing. The modifications to the 
FORCE subroutine are given in Figure B-2, and the subroutine 
CPCALC is listed in Figure B-3. 

The parameters IWOOD and IROW control the calculation of the 
indices of the influence coefficients corresponding to the impact 
theory panels. IROW is the number of panels, leading edge to 
trailing edge, in the linear theory representation. Three options 
are provided for IWOOD. These are: 

IWOOD = 0 The normal operation of the program is unaffected. 

IWOOD = 1 The AIC matrix will be read-in and the combined 
theory pressures, and forces and moments will be 
calculated. The strip input option is used to 
define the geometry. The upper surface panels are 
input in the same order as the linear theory 
panels. The lower surface panels are then input 
in reverse order, i.e., the last panel on the 
lower surface corresponds to the first panel of 
the upper surface. 
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REPLRCE pres S I..IITH 

Cnr-ir-IDH .-' GIiHTR - LTaT !. J » S YMFCT !. I DPN , I GTYPE ? L j I !. I RDI.I 


REPLRCE PRES 53 WITH 

RERD (TRPE I N ? £ U NCnMP, I FSRVE ■■ T I TLE j I WODD < I ROW 

REPLRCE PRES 54 WITH 
£0 FORr-IRTaS, 1173X7 15R4 7£I£> 


Figure B-1. Modifications to PRES Subroutine. 


REPLRCE FDRC 1£ WITH 

CDMHnH •■• GriRTR. -LTnT 7 J 7 SYMFCT 7 I DPh 7 I GTYPE 7 L 7 I WDDD 7 I ROW 


IHSERT BETWEEN FDRC £3 RND FORC £4 

COr'INDN ■••■CPRRRY - ICPS7LTS7CPS<££4;:' 7 BETRM 7 BELT •::££4;:' 

INSERT BETWEEN FORC 35 RND FORC 36 
BETRM=SQPT >:;MRCH**£-1 . > 

ICPS=-1 

1 ICPS=ICPS+1 


INSERT BETWEEN FDRC ££3 RND FDRC ££4 
BELT CLi' =DELTRR 


INSERT 

BETWEEN FDRC 

571 RND FDRC 5 

IF 

acps. EQ. n> 

CPS <0 =CP 

IF 

a CPS. EQ. 1 :• 

cp=cps (;l> 

INSERT 

BETWEEN FDRC: 

749 RND FDRC 7' 

IF 

acPS. EQ. rni 

.RND. aWDDD.NE. 

IF 

aCPS. EQ. l :- 

GD TD 1 


CRLL CPCRLC 


Figure B-2. Modifications to Subroutine FORCE (Next, IG4S), 



on non ooo ooo non non nono 


SUBROUTINE CPCALC 


THIS routine uses the INVERTED MATRIA (A) FROM WOOOWARD 
TO ADJUST THE component CP VALUES CALCULATED IN FORCE SUBH 


COMMON /6OATA/LT0T • JU » S Y mF C I ♦ I URN , 1 (i TYPE * LL t I WOOD tNROW 
COMMON /CPARRY/ I CMS » L TS ♦ CPS t 224 ) , BE T AM , F AC T . D£L T ( 224 ) 
COMMON /TAPE/ TAPE1N,TAPE0 t 

dimension CP. (224) t A (1 10 t UP) t ) »S ( 1 1 0) 

INTEGER TAPEINtTAPEOT 


REAU (A) ARRAY 

NTAPE«iS 
PEWINOIS 
L T=2* (LIS/2) 

L»LT/2 

ho Vo 1=1 »L 

READ (NTAPE,200) ( A ( I » J ) t J» 1 » L » 

10 CONTINUE 

CALCULATE THICKNESS FACTOR 


IS 


DO 15 J*L 1 *LT 

M*NHU«<4 2*NH0w« ( ( J-Ll ) /NRUW) - ( J-Li ) 

IF < I wOOO.ECa, 2 ) Maj-NHUW-2*NROW* t ( j-L 1 ) / NROW) 
£TA« (UELT (M) ♦DELT ( J) ) /2. 

ALPHS (CELT ( J) -OELT (M) ) /2, 

B (M) =l ,♦ 1 ,2*0ETAM*ETA+ .boBETAM**! ,5*ETA0*2 


CALCULATE UPPER SURFACE 


DO 20 1=1, L 
CP ( I ) *0. 

DO 20 J»1,L 

CP<I)*CP(1) A ( I , J) *CPS ( J) 
20 CONTINUE 

CALCULATE LOwEr SURFACE 


DO 30 1=L1,LT 
CP(I)=0, 

K«NROW+2»NROW* ( ( I-Ll ) /NROw) ♦!- ( 1-L ) 

IF ( IwOOO.EQ.2) (V = I-NR0W-2 *nR0w* ( ( 1-L 1 ) /NROW) 

00 25 J=L1»LT 

mcNR0W*2«nr0w* ( ( J-L l ) /NR0 w)-( J-L l ) 

IF ( I WOOO.EU. 2 ) M« J-NR0W-2«nR0w« ( ( j-L 1 ) /NROW ) 

CP(I)=CP(1) * A(K,M)«CPS(J) 

25 CONTINUE 

CP ( I ) =CP ( 1 ) *aErAM/A, 

30 CONTINUE 

WRITE CPS array from FORCE SU0R, 
write (TAPEOT,210) 

WR I TE (TAPEUT ,200) (CPS ( X) , 1= I ,LlS) 

00 40 1=1, LT 
CPS ( 1 ) =CP ( I ) 

40 CONTINUE 

write CP array CALCULATED HERE 
WRITE (lAPeoT,220) 

WRITE (TaPEOT, 200) ( CPS ( I) , 1 = I , L t S ) 

ICPS=1 

WRITE ( TAPEOT,23o) 

wR I TE (T APEOT ,200 ) (r>(I),l = l,L) 

200 FOftflAT ( lX,aElb,8) 

210 format ( IHI ,20X , 3SHC0MP0NENT CP VALUES FROM FORCE SU0R,//) 
220 format ( lHl,2oA,28HADjUSieO COMPONENT CP VALUES,//) 

230 FORMAT(lHl) 

RETURN 

END 


Figure B-3. Subroutine CPCALC 


IWOOD = 2 The combined theory calculations are performed, 
but the geometry is input in the NASA wave drag 
format. 

When the combined calculations are performed, the unmodified 
impact pressures are computed, along with the forces and moments, 
and output in the normal output format. Then the combined theory 
pressures, and the corresponding forces and moments, are calcu- 
lated and output at the same angle of attack. The output is, 
again, in the normal output format. 
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